/* USER CODE BEGIN Header */
/**
  ******************************************************************************
  * @file           : main.c
  * @brief          : Main program body
  ******************************************************************************
  * @attention
  *
  * Copyright (c) 2025 STMicroelectronics.
  * All rights reserved.
  *
  * This software is licensed under terms that can be found in the LICENSE file
  * in the root directory of this software component.
  * If no LICENSE file comes with this software, it is provided AS-IS.
  *
  ******************************************************************************
  */
/* USER CODE END Header */
/* Includes ------------------------------------------------------------------*/
#include "main.h"
#include "adc.h"
#include "dma.h"
#include "tim.h"
#include "usart.h"
#include "gpio.h"

/* Private includes ----------------------------------------------------------*/
/* USER CODE BEGIN Includes */
#include "arm_math.h"
#include "fdacoefs.h"
#include "tmwtypes.h"
#include "stdlib.h"
/* USER CODE END Includes */

/* Private typedef -----------------------------------------------------------*/
/* USER CODE BEGIN PTD */
typedef struct _Biquad{
	real64_T b0, b1, b2;
	real64_T a1, a2;
	real64_T x1, x2;
	real64_T y1, y2;
}Biquad;
typedef struct _IIRFilter{
	Biquad sections[MWSPT_NSEC];
}IIRFilter;
/* USER CODE END PTD */

/* Private define ------------------------------------------------------------*/
/* USER CODE BEGIN PD */

/* USER CODE END PD */

/* Private macro -------------------------------------------------------------*/
/* USER CODE BEGIN PM */

/* USER CODE END PM */

/* Private variables ---------------------------------------------------------*/

/* USER CODE BEGIN PV */
volatile uint32_t gADCValue = 0;
volatile uint8_t gADCState = 0;
const int NL[MWSPT_NSEC] = { 1,3,1,3,1,3,1,3,1 };
const real64_T NUM[MWSPT_NSEC][3] = {
  {
    0.01475498578667,                 0,                 0 
  },
  {
                   1,                 2,                 1 
  },
  {
    0.01360061927639,                 0,                 0 
  },
  {
                   1,                 2,                 1 
  },
  {
    0.01283223764635,                 0,                 0 
  },
  {
                   1,                 2,                 1 
  },
  {
    0.01245152653961,                 0,                 0 
  },
  {
                   1,                 2,                 1 
  },
  {
                   1,                 0,                 0 
  }
};
const int DL[MWSPT_NSEC] = { 1,3,1,3,1,3,1,3,1 };
const real64_T DEN[MWSPT_NSEC][3] = {
  {
                   1,                 0,                 0 
  },
  {
                   1,    -1.84911090872,   0.9081308518666 
  },
  {
                   1,                 0,                 0 
  },
  {
                   1,   -1.704444438844,   0.7588469159497 
  },
  {
                   1,                 0,                 0 
  },
  {
                   1,   -1.608150015067,   0.6594789656524 
  },
  {
                   1,                 0,                 0 
  },
  {
                   1,   -1.560438884015,   0.6102449901735 
  },
  {
                   1,                 0,                 0 
  }
};
/* USER CODE END PV */

/* Private function prototypes -----------------------------------------------*/
void SystemClock_Config(void);
/* USER CODE BEGIN PFP */

/* USER CODE END PFP */

/* Private user code ---------------------------------------------------------*/
/* USER CODE BEGIN 0 */
void HAL_ADC_ConvCpltCallback(ADC_HandleTypeDef* hadc){
	if(hadc->Instance == ADC1){
		gADCValue = HAL_ADC_GetValue(hadc);
//		gADCValue = hadc->Instance->DR;
		gADCState = 1;
		HAL_GPIO_TogglePin(LED_GPIO_Port, LED_Pin);
	}
}

void IIRFilter_Init(IIRFilter *filt, const real64_T NUM[MWSPT_NSEC][3], const real64_T DEN[MWSPT_NSEC][3]){
    for (int i = 0; i < MWSPT_NSEC; i++) {
        filt->sections[i].b0 = NUM[i][0];
        filt->sections[i].b1 = NUM[i][1];
        filt->sections[i].b2 = NUM[i][2];
        filt->sections[i].a1 = -DEN[i][1];
        filt->sections[i].a2 = -DEN[i][2];
        filt->sections[i].x1 = filt->sections[i].x2 = 0;
        filt->sections[i].y1 = filt->sections[i].y2 = 0;
    }
}

real64_T  IIRFilter_Process(IIRFilter *filt, real64_T input){
    real64_T output = input;

    for (int i = 0; i < MWSPT_NSEC; i++) {
        Biquad *bq = &filt->sections[i];

        real64_T yn = bq->b0 * output + 
                      bq->b1 * bq->x1 + 
                      bq->b2 * bq->x2 + 
                      bq->a1 * bq->y1 + 
                      bq->a2 * bq->y2;
				bq->x2 = bq->x1;     
				bq->x1 = output;     
				bq->y2 = bq->y1;     
				bq->y1 = yn;          
        output = yn;
    }

    return output;
}
/* USER CODE END 0 */

/**
  * @brief  The application entry point.
  * @retval int
  */
int main(void)
{
  /* USER CODE BEGIN 1 */
	volatile int ADCVref = 10;
	static float ADCVoltage = 0;
	static float FFTInputData[256*2] = {0};
	static float FFTOutputData[256] = {0};
	arm_cfft_radix4_instance_f32 scfft;
	
	int FrequencyMeterLine[4];
	FrequencyMeterLine[0] = 0x2A;
	FrequencyMeterLine[3] = 0x2B;
	
	int * DynamicArray = (int *)malloc(6 * sizeof(int));
	if(DynamicArray == NULL){
		Error_Handler(); 
	}
	DynamicArray[0] = 10;
	DynamicArray[1] = 20;
	DynamicArray[2] = 30;
	DynamicArray[3] = 40;
	DynamicArray[4] = 50;
	DynamicArray[5] = 60;
	DynamicArray[6] = 70;
  /* USER CODE END 1 */

  /* MCU Configuration--------------------------------------------------------*/

  /* Reset of all peripherals, Initializes the Flash interface and the Systick. */
  HAL_Init();

  /* USER CODE BEGIN Init */

  /* USER CODE END Init */

  /* Configure the system clock */
  SystemClock_Config();

  /* USER CODE BEGIN SysInit */

  /* USER CODE END SysInit */

  /* Initialize all configured peripherals */
  MX_GPIO_Init();
  MX_DMA_Init();
  MX_ADC1_Init();
  MX_TIM3_Init();
  MX_USART1_UART_Init();
  /* USER CODE BEGIN 2 */
	HAL_ADC_Start_IT(&hadc1);
	HAL_TIM_Base_Start(&htim3);
	
	IIRFilter iir;
	IIRFilter_Init(&iir, NUM, DEN);
  /* USER CODE END 2 */

  /* Infinite loop */
  /* USER CODE BEGIN WHILE */
  while (1)
  {
    /* USER CODE END WHILE */

    /* USER CODE BEGIN 3 */
//		for(int i = 0; i < 256; i++){
//			if(gADCState){
//				gADCState = 0; 
//				ADCVoltage = ((float)gADCValue * ADCVref)/(1 << 12)*1000;//mV
//				ADCVoltage = (float)IIRFilter_Process(&iir, (double)ADCVoltage);
//				FFTInputData[i*2] = ADCVoltage;
//				FFTInputData[i*2+1] = 0;
//	//			HAL_UART_Transmit(&huart1, (uint8_t *)&gADCValue, sizeof(gADCValue), 100);
//			}
//			while(!gADCState);
//		}
//		
//		HAL_ADC_Stop_IT(&hadc1);
//		
//		arm_cfft_radix4_init_f32(&scfft, 256, 0, 1);
//		arm_cfft_radix4_f32(&scfft,FFTInputData);
//		arm_cmplx_mag_f32(FFTInputData,FFTOutputData,256);
//		
//		for(int i=1; i<30; i++){
//			float DataTemp = FFTOutputData[i];
//			FrequencyMeterLine[1] = DataTemp;
//			FrequencyMeterLine[2] = i;
//			for(int j=0; j<100; j++){
//				HAL_UART_Transmit(&huart1, (uint8_t *)FrequencyMeterLine, sizeof(FrequencyMeterLine), 100);
//				//while(huart1.gState != HAL_UART_STATE_READY);
//			}
//		}
//		HAL_ADC_Start_IT(&hadc1);
  }
  /* USER CODE END 3 */
}

/**
  * @brief System Clock Configuration
  * @retval None
  */
void SystemClock_Config(void)
{
  RCC_OscInitTypeDef RCC_OscInitStruct = {0};
  RCC_ClkInitTypeDef RCC_ClkInitStruct = {0};
  RCC_PeriphCLKInitTypeDef PeriphClkInit = {0};

  /** Initializes the RCC Oscillators according to the specified parameters
  * in the RCC_OscInitTypeDef structure.
  */
  RCC_OscInitStruct.OscillatorType = RCC_OSCILLATORTYPE_HSE;
  RCC_OscInitStruct.HSEState = RCC_HSE_ON;
  RCC_OscInitStruct.HSEPredivValue = RCC_HSE_PREDIV_DIV1;
  RCC_OscInitStruct.HSIState = RCC_HSI_ON;
  RCC_OscInitStruct.PLL.PLLState = RCC_PLL_ON;
  RCC_OscInitStruct.PLL.PLLSource = RCC_PLLSOURCE_HSE;
  RCC_OscInitStruct.PLL.PLLMUL = RCC_PLL_MUL9;
  if (HAL_RCC_OscConfig(&RCC_OscInitStruct) != HAL_OK)
  {
    Error_Handler();
  }

  /** Initializes the CPU, AHB and APB buses clocks
  */
  RCC_ClkInitStruct.ClockType = RCC_CLOCKTYPE_HCLK|RCC_CLOCKTYPE_SYSCLK
                              |RCC_CLOCKTYPE_PCLK1|RCC_CLOCKTYPE_PCLK2;
  RCC_ClkInitStruct.SYSCLKSource = RCC_SYSCLKSOURCE_PLLCLK;
  RCC_ClkInitStruct.AHBCLKDivider = RCC_SYSCLK_DIV1;
  RCC_ClkInitStruct.APB1CLKDivider = RCC_HCLK_DIV2;
  RCC_ClkInitStruct.APB2CLKDivider = RCC_HCLK_DIV1;

  if (HAL_RCC_ClockConfig(&RCC_ClkInitStruct, FLASH_LATENCY_2) != HAL_OK)
  {
    Error_Handler();
  }
  PeriphClkInit.PeriphClockSelection = RCC_PERIPHCLK_ADC;
  PeriphClkInit.AdcClockSelection = RCC_ADCPCLK2_DIV6;
  if (HAL_RCCEx_PeriphCLKConfig(&PeriphClkInit) != HAL_OK)
  {
    Error_Handler();
  }
}

/* USER CODE BEGIN 4 */

/* USER CODE END 4 */

/**
  * @brief  This function is executed in case of error occurrence.
  * @retval None
  */
void Error_Handler(void)
{
  /* USER CODE BEGIN Error_Handler_Debug */
  /* User can add his own implementation to report the HAL error return state */
  __disable_irq();
  while (1)
  {
  }
  /* USER CODE END Error_Handler_Debug */
}

#ifdef  USE_FULL_ASSERT
/**
  * @brief  Reports the name of the source file and the source line number
  *         where the assert_param error has occurred.
  * @param  file: pointer to the source file name
  * @param  line: assert_param error line source number
  * @retval None
  */
void assert_failed(uint8_t *file, uint32_t line)
{
  /* USER CODE BEGIN 6 */
  /* User can add his own implementation to report the file name and line number,
     ex: printf("Wrong parameters value: file %s on line %d\r\n", file, line) */
  /* USER CODE END 6 */
}
#endif /* USE_FULL_ASSERT */
